Thermo-mechanical behavior of surface acoustic waves in ordered arrays of nanodisks 
studied by near infrared pump-probe diffraction experiments 
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The ultrafast thermal and mechanical dynamics of a two-dimensional lattice of metallic nano- 
disks has been studied by near infrared pump-probe diffraction measurements, over a temporal 
range spanning from 100 fs to several nanoseconds. The experiments demonstrate that, in these 
systems, a two-dimensional surface acoustic wave (2DSAW), with a wavevector given by the recip- 
rocal periodicity of the array, can be excited by ~120 fs Thsapphire laser pulses. In order to clarify 
the interaction between the nanodisks and the substrate, numerical calculations of the elastic eigen- 
modes and simulations of the thermodynamics of the system are developed through finite-element 
analysis. At this light, we unambiguously show that the observed 2DSAW velocity shift originates 
from the mechanical interaction between the 2DSAWs and the nano-disks, while the correlated 
2DSAW damping is due to the energy radiation into the substrate. 
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I. INTRODUCTION 

In recent years, the progress in the fields of transduc- 
ers and sources of coherent acoustic excitations in the 
GHz range, triggered extensive studies on the thermo- 
mechanical properties of solid-state mesoscale systems. 
In particular, time-resolved reflectivity experiments, per- 
formed on gratings of metallic nanometric stripes (2- 
D confined) on transparent (SiO^) or semitransparent 
(Si) substrates, evidenced oscillations in the GHz range 
EL S S 0] > however their origin of these oscillations is 
still matter of debate 0, 13 . The reason resides in the 
difficulty to model the strong mechanical coupling [H, Q 
between the nano-objects and the substrate and to calcu- 
late the mechanical eigenfrequencies of the coupled sys- 
tem. 

In addition, there is a lack of time-resolved experi- 
ments on the thermo-mechanical properties of 3-D con- 
fined nano-objects, mainly because of the difficulties in 
developing an experimental set-up with the required res- 
olution and sensitivity. At the best of our knowledge, the 
only experiment reported in the literature concerns time- 
resolved reflectivity performed on arrays of gold nano- 
metric disks on a S1O2 substrate where the observed 
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variations of the reflectivity signals were attributed to 
the normal modes of the single metallic nano-objects. 

In a recent letter [7[, focused on the study of the 
magnetic dynamics of ordered arrays of magnetic nano- 
disks on silicon, we suggested the possibility of trigger- 
ing a collective acoustic excitation of the nano-disk ar- 
ray/silicon coupled system. To demonstrate that two- 
dimensional surface acoustic waves (2DSAWs) can be ex- 
cited in ordered arrays of metallic nanodisks, by sub-ps 
light pulses, the infrared diffraction time-resolved exper- 
iment reported in Ref. |7j is here further extended and 
improved. In the present work, by measuring the fre- 
quencies of the oscillations of the first order diffracted 
light, as a function of the array periodicity D, we demon- 
strate that two-dimensional SAWs, with a wavevector q 
defined by |g|=27r/|D|, can be excited by a near-infrared 
ultrashort laser pulse. In addition, the comparison be- 
tween finite-element analysis and accurate measurements 
of the 2DSAW velocity shift and damping, as a function 
of q, clearly indicates that the observed 2DSAW velocity 
shift originates from the thermo-mechanical interaction 
between the 2DSAWs and the nano-disks, while the cor- 
related 2DSAW damping is due to the energy radiation 
into the substrate. 

The manuscript is organized as follows: Sec. HH the de- 
tails of the time-resolved technique are presented. In Sec 
IIII1 we report on the time-resolved diffraction measure- 
ments on samples with periodicities of 2000 nm, 1 000 nm, 
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FIG. 1: A scheme of the time-resolved optical set-up is 
reported. The variation of the intensity of the first-order 
diffracted probe beam is measured through a lock-in amplifier 
referenced at 100 kHz. 

800 nm and 600 nm and the results are compared to the 
calculations of the mechanical eigenmodes. In Sec. IIVI 
we discuss both the numerical simulations of the com- 
plete thermal evolution of the nanodisks/substrate cou- 
pled system and the model developed to account for 
the measured dynamics of the intensity variation of the 
diffracted beam. In Sec. |V] we address the mechanism 
responsible for the 2DSAW velocity shift and damping. 



II. EXPERIMENTAL SET-UP 

The time-resolved optical set-up, sketched in Fig. [TJ 
has been developed to achieve a high-sensitivity in very 
reduced perturbing conditions. 

The light source is a 76 MHz repetition rate 
Ti:Sapphire oscillator, operating at a central wavelength 
of ~790 nm and with a pulse timewidth (At) of 120 fs. 
The 2 W laser beam is divided in the pump and probe 
beams. Both the pump and probe intensities are con- 
trolled by means of a half- wavelength (A/2) plate and a 
calcite polarizer. The energy of the pump is <10 nJ, 
while less than 1 nJ is used for the probe. The pump 
beam is modulated at 100 kHz, through a photo-elastic 
modulator (PEM) placed between two crossed polarizers. 
To optimize the optical signal to noise ratio, the spot 
sizes of the pump and probe beams was set at 60 /im 
and 40 /j,m, respectively. Second harmonic generation 
by a BBO crystal was used for the temporal and spatial 
overlap of the pump and probe beams. 

The first-order diffracted probe beam, induced by the 
periodicity of the system, is focused and detected by a Si- 
photodiode. The noise due to the direct detection of the 
pump beam scattered light is reduced through a crossed- 
polarizer, placed in front of the detector. The detected 
signal is filtered by a lock-in amplifier referenced to the 
pump modulation frequency. The phase of the signal 
is measured by locking the internal phase of the lock-in 



FIG. 2: The detail of the pump beam modulation, through the 
PEM between two crossed polarizers, is shown. A modulation 
of the intensity of the pump beam at 100 kHz is obtained. On 
the right the unit cell of the nanodot array is represented. 

amplifier to the intensity modulation of the pump beam. 
In this way it is possible to acquire signals in-phase and 
7r/2-out-of-phase with respect to the pump modulation. 

To investigate the relaxation dynamics in the nanosec- 
ond time-domain, a long range translator is used to 
change the pump optical path up to 1.2 m, corresponding 
to a maximum delay between the pump and probe pulses 
of ^4 ns. Two spatial points of the pump beam are im- 
aged onto two quadrant photodiodes interfaced to a feed- 
back system that drives two piezo-nano motors mounted 
on optical mirrors. This system is used to keep the spatial 
position of the pump beam fixed during the experiment 
and avoid variations of the detected signal due to the 
pump-probe misalignment. 

Permalloy (Py) nanosized cylinders (nanodisks), ar- 
ranged in an ordered square lattice, have been syn- 
thetized by electron-beam lithography and lift-off tech- 
niques on a Si(100) surface. The array periods, the nan- 
odisk thicknesses and diameters have been carefully mea- 
sured by AFM (atomic force microscopy) and grazing- 
angle X-ray reflectivity. Time-resolved opical measure- 
ments have been performed on nano-disk arrays with 
a nominal thickness of 30 nm and periods of 600 nm, 
800 nm, 1000 nm and 2000 nm. The diameter of the 
disks, as measured by AFM, is half of the array period. 
A Peltier device is used to keep constant the sample back- 
side temperature during the experiment. 

III. TIME-RESOLVED REFLECTIVITY 
MEASUREMENTS 

In Fig. [31 the relative variation (AI\d/Iid) of the in- 
tensity of the beam diffracted from the 800 nm-period 
sample is reported in the sub-nanosecond timescale. 
To decouple the temperature variation in the sub- 
nanosecond timescale from the average heating of the 
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FIG. 3: Transient reflectivity on permalloy nanodisks of 
800 nm period and 400 nm diameter at different pump flu- 
ences. In the inset the dependence of the background level 
and of the oscillation amplitude is reported as a function of 
the pump intensity. 



substrate, only the in-phase component of the variation 
is displayed, see Appendix |A"1 I\u is the measured total 
diffracted intensity. Several sets of measurements were 
performed by changing the intensity of the pump excita- 
tion. In particular, we varied the average incident power 
between 225 mW (energy: 2.9 nJ/pulse) and 640 mW 
(energy: 8.4 nJ/pulse). 

After pump excitation at t—0, the signal quickly rises 
and starts to oscillate with a frequency of ^175 ps. The 
variation of the diffracted intensity is of the order of 
A7i£)//i£)~l-5T0 _5 and positive, i.e. an increase of the 
intensity is detected after excitation. At negative delays, 
i.e. when the probe beam investigates the sample before 
excitation by the pump pulse, a positive background is 
detected. This signal is attributed to the mean heating 
of the dot/silicon system related to the high repetition 
rate of the laser beam, as discussed in detail in Sec. IIVI 

In the inset of Fig. O the background level (triangles) 
and the maximum oscillation amplitude (circles) are re- 
ported as a function of the incident power. The solid lines 
are the linear fit to data. These results indicate that both 
the background and the oscillation linearly depend on the 
excitation energy and, as a consequence, on the temper- 
ature increase of the metallic array. This demonstrates 
that the temperature increase is small enough to prevent 
non-linearities caused by the temperature variation of the 
relevant optical, thermal and mechanical constants of the 
system. 

To investigate the physical origin of the measured os- 
cillation, we extended the experiment to the nanosecond 
timescale and measured samples with different periodic- 
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FIG. 4: Long-range transient reflectivity (colored circles) on 
permalloy nanodisks arrays with different periodicity are re- 
ported. The lattice period D and disk diameter 2a are re- 
ported for each sample. The data has been fitted (solid lines) 
with Eq. [H 



ities and the same disk diameter to period ratio. From 
the results reported in Fig. [4] we note that both the 
oscillation frequency and damping change as the peri- 
ods are varied. In particular, we measured a variation 
of the oscillation frequencies from ~409 ps to ^135 ps, 
by changing the nominal periodicity from 2000 nm to 
600 nm. In Fig. [5] we display the oscillation frequencies 
of the measurements reported in Fig. [4] (yellow circles), 
as a function of the array periodicity. Assuming that 
the wavelength of the displacement is given by the pe- 
riodicity of the system, the wave velocity vs is obtained 
performing a fit of the the data to a power function with 
the exponent n=-l. The result is v$ — 4800 ± 100 m/s. 
This value is ~ 5% below the reported SAW velocities 
on [001] Si surfaces (4921 m/s in the (100)direction and 
5070 m/s in the (110) direction, see Ref. Q). Moreover, 
the small, systematic shift in the measured velocity with 
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respect to the expected SAW velocity of a pure [001] Si 
surface can be explained by the mechanical-loading due 
to the permalloy disks, as will be shown below. These 
facts demonstrate that the oscillations in the data of Fig. 
|4] are caused by 2DSAWs with different q on the Si sur- 
face. 

To fully support this interpretation, the eigenmodes 
of the single (decoupled) nanodisks have to be excluded 
as a possible origin for the oscillations. To do so, we 
performed a finite element analysis of the lowest three 
eigenfrequencies of single nanodisks and of the 2DSAW 
modes, taking into account the mechanical loading of the 
nanodisks. The measured sample dimensions have been 
used throughout the simulations. 

The calculated eigenmodes are presented in Figj5] The 
shape reproduces, linearly expanded, the real shape of 
the modes, whereas the color scale, decreasing from red 
to blue, refers to the modulus of the total displace- 
ment. Note that these calculations do not require any ad- 
justable parameters (the permalloy stiffness coefficients 
have been calculated from the Ni and Fe polycrystalline 
coefficients given in the Appendix of Ref. M). In the 
case of the single nanodisk eigenfrequencies, both the val- 
ues of the frequencies and the behavior as a function of 
the disks diameter differ markedly from the experimen- 
tal data. On the contrary, the simulation of the 2DSAW 
reproduces within few percents the experimental frequen- 
cies. We can therefore unambiguously conclude that the 
observed oscillations are caused by 2DSAW propagating 
on the Si surface. 

At the light of the results reported in this section, 
we can rationalize the physics of the system as fol- 
lows. The laser-induced impulsive heating is at the ori- 
gin of the disks expansion. The fast expansion of the 
metallic nanostructures, mechanically strongly coupled 
to the substrate, launches the 2DSAW in the silicon. 
On the subnanosecond timescale, the 2DSAW modulates 
the disk radius Sa(t), inducing the measured oscillation 
of the diffracted signal Aim. In the following section 
we develop a model to quantitatively interpret the time- 
dependence of Aii£>. 



IV. THERMODYNAMIC AL AND 
MECHANICAL MODELS 

A. Thermodynamical model 

The absorption of a sub-ps laser pulse induces an im- 
pulsive heating of the system. The larger increase of the 
temperature of the metallic nanodisks with respect to 
the substrate, triggers a spatially modulated stress on 
the silicon surface. The mechanical problem under inves- 
tigation is therefore entangled with the thermodynamics 
of the system. The detailed understanding of thermo- 
mechanical dynamics involved in the experiments, relies 
on the knowledge of the spatial and temporal tempera- 
ture profiles. 
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FIG. 5: Comparison of the measured oscillation frequencies 
(yellow circles) to the calculated eigenfrequencies (black tri- 
angles) of the 2DSAW in silicon coupled with the metallic 
nanodisks, as a function of the array periodicity. The solid 
yellow line is the best fit of the data to a power function. The 
determined power exponent is n=-l. In addition, Fig. [S] dis- 
plays the calculated lower three eigenfrequencies (blue dots) 
of the isolated disks, that is decoupled from the substrate, as 
a function of the nanodisks radius. In all the reported cal- 
culation the measured thicknesses of the samples have been 
taken into account. 



A numerical simulation of the excitation and heat- 
diffusion process, shows that three related mechanisms 
must be considered. First, the average heating occurs on 
the milliseconds time scale and it can therefore be exper- 
imentally filtered out. Second, a nano-disk temperature 
increase of ~ 10 K is attained on the picosecond time- 
scale and the related thermal expansion is responsible for 
launching the 2DSAW. Third, the disk-substrate heat ex- 
change process is effective on the nanosecond time-scale. 

The numerical values of the physical constants used 
in this section are reported in Table IIII1 Appendix [C] 
Note that, in the present calculations, the boundaries 
of the disks are not allowed to expand, which limits the 
results to the thermal problem. A full thermo-mechanical 
analysis, which is out of the scope of the present paper, 
requires the inclusion in the simulations of the thermal 
expansion term. 

Modeling the average heating requires taking the en- 
tire sample volume into account. The equilibrium tem- 
perature distribution, and the timescale to attain it, are 
found solving the thermal problem for a continuous laser 
(switch-on rise-time ~100 fs, average power ~0.4 W and 
beam spot radius ~27.5 fim) shone on a silicon cylinder 
of the size of the substrate used in the actual experiment 
(5 mm radius x 1 mm height; the dots volume can be 
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FIG. 6: Equilibrium temperature profile of the Si slab at- 
tained 1 ms after the arrival of the first laser pulse. The tem- 
perature profile is shown from the surface down to a depth of 
100 j-im, vertical axis, and within a 100 fim radius from the 
Si substrate center, horizontal axis. 

neglected, being 2-1CP 6 of the total volume): 

C Sl (T Si ) = Peon + k Si V 2 T Sl (1) 

where Tsi, Csi, k$i are the Si lattice temperature, 
heat capacity and diffusion coefficient respectively. The 
power per unit volume absorbed at time t and depth 
1 2 1 by the Si sample from the continuous laser beam is 
P C on(t,z,r)=P con (t)P con (z,r) (see appendix [C] for the de- 
tailed analytical expression). 

As for the boundary conditions, the bottom surface 
of the Si cylinder is set at the constant temperature of 
293.5 K, i.e. the temperature of the Peltier cell holding 
the sample. The heat flux boundary condition is imposed 
on all other surfaces: 

n St ■ kstfTsi = h(T air - T Sl ) (2) 

nsi being the unit vector normal to the substrate surface, 
T a i r — 293.5 K the ambient air temperature, h the Si-air 
heat transfer coefficient. The equilibrium temperature 
distribution is established after ~1 ms and is shown in 
Fig. O The portion of the Si slab exposed to the laser 
beam has a rather constant temperature down to a depth 
of 10 fiTa, the average temperature in this volume being 
-305 K. 

The obtained result is useful for two main reasons. 
First, the time the period of the PEM modulation 
(T=l/100 kHz=10 /us) of the probe is at least 100 times 
smaller than the time necessary to heat the entire sam- 
ple volume. As a consequence, we can filter out the 
signal related to the average heating by acquiring only 
the component of the diffracted signal in-phase with the 
pump modulation (see Appendix |A"]) . Second, in further 
determining the picosecond dynamical evolution of the 
electron and lattice temperatures of the system, we will 
assume an initial dot and Si temperatures of 305 K. 

After excitation by a single pump pulse, the thermal 
evolution problem spans three time scales. In the first 
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FIG. 7: The simulation of the time evolution of the average 
electron temperature (T e ) and average lattice temperature of 
the permalloy disk (Tp y ) is shown. The maximum tempera- 
ture of the electron system (~350 K) is reached 150 fs after 
laser excitation, whereas the thermalization with the lattice is 
completed on the picosecond timescale. An equilibrium tem- 
perature of 313 K is established ~2 ps after sample excitation. 
The left inset displays the permalloy lattice temperature time 
evolution at different heights within the dot. The heights are 
expressed in nm with respect to the Si surface. In the right 
inset the time-resolved reflectivity measurements, performed 
on thin film samples, are reported. 



step, the laser short pulse heats the electron gas of the 
metallic nanodisks (sub-picosecond time scale). In the 
second step, the hot electron gas thermalizes with the lat- 
tice (picosecond timescale). In the third step, the disks 
thermalize with the silicon substrate (nanosecond time 
scale). This three-steps sequence repeats itself upon ar- 
rival of a new laser pulse each 13 nsec. 

Finite-element analysis was performed on the unit cell 
reported in Fig. [2j This simplification is justified a- 
posteriori, the simulations showing that there is no ther- 
mal interaction among adjacent cells, i.e. the tempera- 
ture asymptotic values are established well within the cell 
boundaries. The unit cell boundary is actually taken as 
cylindrical, for sake of simplicity in writing the equations. 

The time evolution of the electron and lattice temper- 
atures in the nanodisks is given by the coupled equations 
system [loj : 

dT 

Ce(T e )-^ = -G • (T e - T Py ) +P p + k e v 2 T e (3) 

Cpy(Tp y ) — = G ■ (T e — Tpy) + kp y XJ Tpy 

where T, k and C indicate temperature, thermal diffusion 
coefficient and heat capacity, respectively, the reference 
to the electron (e) or Permalloy (Py) being indicated by 
the subscript. G is the electron-phonon coupling constant 
and C e = "fT e . The profile of the pulsed power per unit 
volume absorbed by the sample is P p (t,z,r)=P p (t)P p (z,r). 
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The analytical expression for P p (t,z,r) is reported in Ap- 
pendix [UJ 

As for the boundary conditions, the dot is thermally 
coupled to the underlying Si substrate via a thermal 
boundary resistance. The simulations show that the ther- 
mal coupling timescale is ~1 ns, whereas the electron- 
phonon thermalization within the dot occurs on the pi- 
cosecond time-scale. For this reason we discuss the 
sub-ps laser-electrons interaction and the ps electrons- 
phonons thermodynamics on the basis of equations ([3|) 
only, assuming the boundary with Si to be thermally in- 
sulated. The dot free boundaries are coupled to the sur- 
rounding atmosphere via convection, that is the bound- 
ary condition of the kind reported in Eq. [2] applies, after 
the obvious substitutions for the dot electron and lattice 
temperatures, thermal diffusion, and unit vector defining 
the dot boundaries. 

Simulations results are reported in Fig. [7J The laser 
pulse is absorbed by the electrons, that reach a maxi- 
mum temperature of ^350 K in the first few hundreds 
of femtoseconds. The electron gas then thermalizes with 
the permalloy lattice ~2 ps after laser excitation, the 
two populations attaining an asymptotic temperature 
^313 K. The permalloy lattice temperature evolution is 
calculated at different heights within the dot, the heights 
being expressed with respect to the Si surface. The max- 
imum lattice temperature spread is ~4.5 K and occurs 
during the electron-lattice thermalization process, see the 
left inset of Fig. [7] This small spread is related to 
the finite thermal diffusion of the electrons in permalloy, 
smaller than in common metals. 

To test the reliability of the thermodynamic simula- 
tions, we performed time-resolved reflectivity measure- 
ments on permalloy films of different thicknesses, de- 
posited on a silicon substrate. In this case, due to the 
lack of spatial modulation, no SAWs are induced and the 
electron-phonon relaxation dynamics is more easily ac- 
cessible. The results are reported in the right inset of 
Fig. El for permalloy thicknesses of 25 nm and 12.5 nm. 
These data can be well described by Eq. ((3]), assuming 
an electron-phonon coupling constant G of the order of 
10 18 W/ (m 3 K). This value is larger than the values mea- 
sured on standard metals, such as gold and silver, but in 
agreement with the results previously reported for pure 
nickel [ll|. It is important to note that, while the reflec- 
tivity variation at negative delays and at delays larger 
than 2 ps is linear with the film thickness, the same sig- 
nal variation is detected in the first hundreds of femtosec- 
onds, independent of the sample thickness. This is due 
to the fact that the reflectivity signal on the femtosecond 
time-scale depends only on the light penetration depth 
(~15 nm), whereas the slow component in the transient 
signal, attributed to the lattice heating, directly depends 
on the absorbed heat per unit volume, which is inversely 
proportional to the sample thickness. 

At the light of these results, we can conclude that, after 
the excitation by the pump pulse, the temperature of the 
nanodisks is impulsively increased of about 8 K with re- 



spect to the equilibrium temperature. Considering that 
the thermal expansion of the permalloy disk induces a 
displacement of the silicon substrate, we can estimate 
that, in the first picoseconds, the displacement of the disk 
radius coupled to the substrate is 5a/a=a ■ AT, where 
a is the thermal expansion coefficient of the permalloy 
disk/silicon coupled system. Assuming ac^asi, being 
asi=2.Q ■ 10~ 6 K _1 [12J the Si thermal expansion co- 
efficient at room temperature, and the temperature in- 
crease of the metallic nano-object AT ~8 K, we obtain 
(5a/a~2.5-lCT 5 . 

On longer timescales, we are going to consider the 
metallic nanodisks as a system in thermodynamical equi- 
librium, which is slowly thermalizing with the cold sub- 
strate acting as a heat reservoir. In this third step the 
coupling between the isothermal dot and the Si substrate 
comes into play, via the boundary conditions: 

n Py ■ kpyWTpy + (T > y ~ T Si )/R = (4) 

-n Sl ■ ksiVTsi - (T Py - T Si )/R = 

Hpj,, risi being the outward unit vector normal to the 
dot and Si boundary respectively and R the interface 
thermal resistance. The interface thermal resistance is a 
thermodynamical parameter accounting for microscopic 
processes such as phonon dispersion mismatch between 
the disk and Si or the presence of a thin SiC>2 layer be- 
tween the former two materials. The value for the in- 
terface resistance was taken ~ 10 -8 Km 2 W _1 a value 
compatible with values found in the literature for similar 
interfaces [HG3. 

The outcome, relevant to the present discussion, is that 
the thermalization occurs on the ns time scale. The detail 
of the coupling mechanism, however interesting, will be 
discussed elsewhere. 



B. Oscillation model 

In analogy with the displacive excitation of coherent 
phonons [15j, we consider the equation governing the 
time dependence of the dot radius, i.e. a simple harmonic 
oscillator equation, with a varying equilibrium position: 

5a(t) = -u%[5a(t) - 5a a (t)} ~ 2 1 8a{t) (5) 

where loq is the oscillation frequency of the surface wave, 
8ao(t) is the radial variation of a non-oscillating disk, 
depending on the lattice temperature and, as a conse- 
quence, on the delay-time from the excitation pulse, and 
7 is the damping constant of the mode. We assume a di- 
rect proportionality between the equilibrium radial vari- 
ation and the temperature increase: 

5a (*) = la5T{t) (6) 

where I is a length parameter and a is the thermal ex- 
pansion coefficient of permalloy. Starting from the simu- 
lations of the thermal problem described in Sec. IIV1 on 
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the nanosecond timescale we can define a unique tem- 
perature T(t) of the nanodisks and assume the following 
dependence: 



disks and the substrate. Substituting eq. [7] in eq. we 
obtain, in the frequency domain: 



ST(t) = 6T(0)e~ t/T H(t) 



(7) 



8a(u)) 



la8T ui 2 



where ST(0) is the temperature difference between the 
disks and the substrate at zero delay and r is the time 
constant regulating the heat exchange between the nano- 



J 



(1/t + ilu) (uJq - ui 2 + 2ijuj) 
By anti-transforming eq. [HI we obtain: 



(8) 



Ht) = laST^4 [ + °° e -t/T H{t)( 



-i(t-t') 



and, by solving: 



smw(t - t')H(t - t')d£ 



(9) 



Sa(t) = 



laSToojQ 



j 2 + l/r 2 -2 1 /r l 



-t/T 
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e" 7 * cosujt + ^e -7 * sinZJt]H(t) 



(10) 



r 



where w = \J uj^ — 7 2 is the renormalized frequency and 
0=l/r-r 

In Fig. [5J Eq. [TO]is displayed in the limit of very short 
(r=10 ps, red line) and very long (t=10 ns, blue line) 
relaxation time constants. In the first case, 1/t ~ oo 
and a/U 1 and, as a consequence, the time evolu- 
tion of the oscillating mode is dominated by the sin-like 
term, in contrast with the experimental results. For this 
reason we can exclude, as the source of the oscillating 
waves, a role played by the photo-excited electrons, be- 
cause the thermalization of the non-equilibrium electron 
population acts on the femtosecond timescale. On the 
contrary, in the case of very slow relaxation dynamics, 
the cos-like term dominates and a damped oscillation su- 
perimposed to a step function is expected (see Fig. [8j. 
The results reported in Fig. [4] indicate that the 2DSAWs 
are launched by a variation of the disks radius which re- 
laxes on a timescale longer than the oscillation period 
(100-400 ps) but comparable to the timescale probed by 
our set-up. 

On the basis of the previous considerations, we are now 
able to evaluate the intensity variation of the first-order 
diffracted beam, related to the oscillation and relaxation 
of the disk radius: 



AJidW 1 .dim ( .^9I 1D dl u 

: = 7 — (~o — Sn(t) + ——8T{t) + -—— 

id 1id on ill aa 



I 



■Sa(t)) 
(11) 

where Sn(t) is the equilibrium electron distribution vari- 
ation and Sa(t) is given by Eq. [TU] The first term is the 
fast response related to the laser-induced non-equilibrium 



TABLE I: The oscillation periods 2tt/uj, the inverse damp- 
ing I/7, obtained from the fit of Eq. (|12|1 to data, are re- 
ported. The periods of the arrays are obtained through two- 
dimensional FFT of the AFM images of the samples. 



Periods (nm) 



2n/u (ps) 



1/7 (Ps) 



a) 


610 ±3 


134.8 ± 0.1 


950 ± 30 


b) 


810 ±10 


175 ±0.1 


1690 ± 60 


c) 


1020 ± 50 


211.2 ± 0.1 


3980 ± 300 


d) 


2018 ± 30 


409.6 ± 0.3 


14000 ± 5000 



electron population (see Sec. IIVP and can be neglected. 
The second term is related to the dependence of the op- 
tical constants of the system on the lattice temperature 
and the third term takes into account the oscillation of 
the dot radius. By using Eq. (JB6J) and © in Eq. (fTTj) . 
we obtain: 



= KST(0)e-^H(t) 

I ID 



+ —5a{t) (12) 



being K=(Iiu)~ 1 (dIiD/dT) a coefficient related to the 
dependence of the permalloy reflectivity on the temper- 
ature. Considering the impulsive expansion of the disk 
radius coupled to the substrate, i.e. Ja/a~2.5T0 -5 (see 
Sec. II Vp . we can use Eq. [T2]to estimate the oscillation in- 
duced in the diffracted signal by the 2DSAW excitation. 
The estimated oscillation amplitude 8I\d/ Tid— 6T0 -5 
is compatible with the oscillation detected in the time- 
resolved measurements. 



T=10 PS 




delay(ps) 

FIG. 8: In the figure the fit (green solid line) of Eq. \TU\to the 
time-resolved diffracted intensity variation (green circles) of 
the 300nm-diameter sample is displayed. The red (top) and 
blue (bottom) solid lines are the expected diffracted intensity 
variations for impulsive and step relaxation dynamics, respec- 
tively. The dotted lines represents the temporal evolution of 
the equilibrium radius variation of the metallic disks. 

By fitting the data with the expression given by Eq. 
(fT2")l it is possible to determine two important parameters 
related to the mechanical properties of the nano-object 
array coupled with the substrate: 

• m~ojQ is the oscillation frequency, mainly related 
to the elastic constants of the substrate 

• 7 is the damping of the 2DSAW 

The data reported in Fig. 0] have been fitted with the 
function reported in Eq. (| 1 2[) and the oscillation fre- 
quency and the damping have been determined with an 
uncertainty drastically smaller compared to the results 
previously reported in the literature on similar systems 
HSU. In Table U we summarize the results obtained. 

V. SAW VELOCITY SHIFT AND DAMPING 

2DSAWs propagate in a thin layer on the surface of the 
substrate, the displacement field vanishing in the sub- 
strate within a depth of the order of the wave length D. 
Therefore, at contrary to bulk modes, the presence of an 
overlay, like the permalloy disks, is expected to modify 
the properties of the surface wave. The resulting varia- 
tion of the wave velocity, in particular, is an aspect of 
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<f -0.3 
-0.4 
-0.5 







1 ^ 2SFTT^ 1 

A measured velocity shift 




calculated velocity shift 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
(h/D)T| x 10" 2 



FIG. 9: The measured 2DSAW velocity shift is reported, as 
a function of the parameter (h/D)ri. The dotted line is the 
velocity shift estimated through the perturbative approach. 



TABLE II: Measured VS AW,Meas. and calculated vsAw,Caic. 
SAW velocities. The calculated velocities have been obtained 
using a perturbative approach reported in Refll7l 





D (nm) 


2a (nm) 


h (nm) 


VMeas. (m/s) 


vcaic. (m/s) 


a) 


610 ±3 


320 ± 10 


60 ±20 


4525 


4502 


b) 


810 ± 10 


380 ± 20 


33 ±5 


4629 


4784 


c) 


1020 ± 50 


470 ± 10 


21 ±2 


4834 


4853 


d) 


2018 ± 30 


990 ± 10 


31 ± 1 


4930 


4865 



practical interest for SAW device applications and has 
been investigated by many authors 0, [ll| [T3] • 

We limit the present discussion to a perturbative ap- 
proach presented for example by Robinson et al. [ItJ • In 
this calculation, the velocity shift is proportional to the 
reflection coefficient r s of the 2DSAW by the overlay, i.e. 
Avsaw /vsaw — r s(h/ D)rj, where h is the thickness of 
the disk and r]=ira 2 / D 2 the filling factor, that is approx- 
imately 7r/16 for the samples of Fig. @] In the case of 
permalloy on Si, it leads to r s = —3.81 [l8|], to be com- 
pared to To = —1.26 for Al on Si and r s = —13.6 for 
Au on Si [4]. These values show that the density con- 
trast between the metallic overlay and the Si substrate is 
the dominant factor determining the velocity shift (mass- 
loading effect). The velocities of the 2DSAW, corrected 
for loading effect of the nano-disks, are presented in the 
Table |TT] and Fig. [5J together with experimental values. 
The calculated velocities, for which, again, no free param- 
eter have been used, are in reasonable agreement with the 
measured ones. The variation of the experimental data 
at small (h/D)rj may be explained by the spread of the 
nanodisks thickness in the samples (see Table HI)) . 

The dependence of the damping coefficient 7 is a cru- 
cial parameter to properly address the physical mecha- 
nisms responsible for the attenuation of the 2DSAW. In 
particular, a q A dependence is expected in the case of 
energy radiation from the SAW to the bulk whereas 
jocq 2 for dissipative damping caused by thermoelastic 
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effect, Akhieser effect Q or viscosity of the electron gas 

Q. 

The origin of the energy radiation of 2DSAW is the 
following: for a free, uniform surface, the 2DSAW is an 
eigenmodc of the system and no radiation is expected. 
When an overlay covers the surface, however, a non-zero 
stress is required at the substrate-overlay interface to 
force the overlay to follow the motion of the substrate 
surface. This stress, more precisely its perpendicular 
component, is allowed to radiate to the bulk. We follow 
here a discussion by Lin et al. [l[ to derive the damping 
expected in that case. The perpendicular stress at the 
substrate-overlay interface is derived from the equation 
of motion, which gives (in the limit of thin overlays): 



&zz = PPyhpyUzQUl 2 



(13) 



where pp y and hp y are the density and height of the 
permalloy film and u z0 is the maximum displacement in 
the perpendicular direction. The corresponding energy 
radiation can be derived by applying a simple relation 
for power flow density of acoustic waves [l9[ , which leads 
to the following radiation rate per unit surface A: 



j/A cx 



Ppyh 



PSiVl,Si 



(14) 



where psi and vi t si are the density and longitudinal veloc- 
ity of the Silicon substrate. To obtain a relation between 
the damping and the frequency, we still have to calculate 
the displacement amplitude. An estimate for this value 
can be obtained by considering that, in our system, the 
temperature increase is kept small enough to stay in the 
linear regime. We can therefore assume that u z q, as well, 
depends linearly on the temperature. Considering the 
same absorbed intensity, the temperature variation is in- 
versely proportional to the mass, i.e. to the thickness of 
the disks. In this case, one obtains the simple relation: 



u z q cx AT cx 1/hpy 



(15) 



where AT is the adiabatic temperature increase of the 
disks. Putting together these relations leads to: 



cx A . 
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(16) 



where A is the 2DSAW wavelength equal to the array 
period. 

In Fig llOl the inverse of the experimental damping 
rates, normalized to the nanodisks surfaces, are reported 
as a function of the frequency of the 2DSAW for the four 
different samples. The data can be satisfactory fitted to 
a power function with an exponent rt=4, in very good 
agreement with the relation reported in Eq. 1161 

Two conclusions can be drawn from this result: first, 
the radiation of the 2DSAW to the bulk modes com- 
pletely dominates the dissipative attenuation, which may 
be explained by the fact that the substrate is made of 
semi-conductor and crystalline material. Second, in our 




1000 

array period (nm) 

FIG. 10: The damping normalized to the nanodisks surface is 
reported as a function of the array periodicity. The damping 
has been obtained by fitting the data reported in Fig. 1101 to 
the function given by (|10|l and (|12[) . The solid line is a power 
function with an exponent n=4 



thermo-mechanical system, the two assumptions used in 
the calculation of the damping, i.e. the adiabaticity of 
the disks at initial times and the linear relation of the 
displacement vs. the temperature, result to be valid. 
This is an important result that demonstrate that such a 
system is suitable for thermodynamics measurements of 
nanoparticles. 



VI. PERSPECTIVES AND CONCLUSIONS 

The present results strongly encourage for the devel- 
opment of a more detailed model of the heat exchange 
process between the metallic disks and the substrate, in 
order to deeply investigate the physics related to the time 
constant r regulating this process (see Section IIVBI) . 
As discussed in this work, the mechanics and thermody- 
namics of the system are intimately related. The laser- 
induced impulsive heating of the metallic nanostructures 
triggers both the heat exchange with the substrate and 
the transfer of mechanical energy through excitation of 
the surface acoustic waves. The study of this physics, 
as the size of the nanodisks is reduced down to tens of 
nanometers and the temperature decreased down to a 
few K, could be of particular relevance in the field of the 
nanocalorimetry. In this regime the thermal population 
of the phonons in silicon is quenched and the thermal 
conductance can exhibit the signature of the quantiza- 
tion of the modes of the nanostructures [2(1 Hl[ . It is im- 
portant to underline that time-resolved experiments have 
already been performed on disordered metallic nanopar- 
ticles in an amorphous dielectric 0, [23j]. However, in 
these systems, the development of a detailed model of 
the laser-induced impulsive heating and the following re- 
laxation process is very difficult, as a consequence of the 
not well-defined thermodynamic boundary conditions at 
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the nanoparticles/dielectric interface. 

At the same time the problem of the elastic coupling 
of the SAWs with the array of nanostructures deposited 
on the substrate is worth of investigation. In fact, the 
interaction of the SAWs with the periodic modulation of 
the surface stress can result in a deformation of the dis- 
persion relation near the borders of the super-Brillouin 
zone and in the opening of a gap in the elastic surface 
modes. Recently it has been demonstrated that the for- 
bidden acoustic gaps, typical of the phononic crystals 
[24| . can be extended to the hypersonic (GHz) range in 
colloidal crystals formed by self-assembly of polystyrene 
nanoparticles [25| . However, in this system, the modes 
can be excited only by thermal motion and no external 
control is possible. On the contrary, two-dimensional ar- 
rays of metallic nanostructures, deposited on a substrate, 
are the best candidates as hypersonic phononic crystals 
where the modes can be coherently excited and controlled 
by light. 

Finally, the study of the optical properties of metallic 
nano-objects has received a growing attention in the last 
years. In particular, when the typical size is reduced to be 
smaller than the light wavelength in the visible spectrum, 
the nano-objects can be represented as classical dipoles, 
i.e. optical antennas, interacting with the electromag- 
netic field. The resonant collective excitation of the free 
electron gas is responsible for unusual and interesting 
electromagnetic properties such as subwavelength plas- 
mon resonances [26| . super-continuum generation due to 
strong field enhancement [271 . l28j | , negative permeability 
[29l | and diffraction from ordered sub-wavelength aper- 
tures in metallic surfaces [30l . |3~J | . A deeper understand- 
ing of the physical processes at this length-scale is the 
first step toward the development of devices in which the 
sub-wavelength optics is coupled to the nano-mechanics 
and thermodynamics. 

In conclusion, we have developed a time-resolved op- 
tical spectroscopy in order to investigate the elastic and 
thermodynamical properties of an array of metallic nano- 
objects deposited on a silicon substrate. Exploiting the 
2-d periodicity of the system, we measured, directly in 
the time-domain, the variation of the diffracted inten- 
sity as a function of the delay from the excitation pulse. 
By changing the periodicity of the array we were able 
to demonstrate that surface acoustic waves are excited 
in the substrate. The measurements of the 2DSAW ve- 
locity shift and damping as a function of the dimensions 
of the system, allowed to discuss the coupling between 

I 



the 2DSAW and the nanostructures and the mechanisms 
responsible for the damping. In addition, numerical cal- 
culations of the eigenmodes of the system and simula- 
tions of the thermalization process, after excitation with 
the pump beam, were performed through finite element 
analysis. Good agreement with the measurements is ob- 
tained. These results demonstrate that metallic nano- 
disks arrays on semitransparent substrates are a partic- 
ularly simple model-system to study the mechanics and 
thermodynamics of nanometric objects coupled to a sub- 
strate. 



APPENDIX A: AVERAGE HEATING 
DECOUPLING 

The pump beam is modulated through a photo-elastic 
modulator (PEM). When the PEM is used between two 
crossed polarizers and the maximum amplitude of the 
mechanical stress corresponds to a retardation of half- 
wavelength of the components along the two axis, the 
modulation of the light intensity is given by: 

I(t) = 7 /2[l - cos(TTsin(ujpt))} (Al) 

where Iq is the incident light intensity and ujp=2itf is the 
PEM fundamental frequency The main feature of I(t) is 
that the modulation is at twice the PEM frequency, i.e. 
100 kHz. 

The pump high-frequency modulation is an important 
aspect for the experiment, since the detector 1// noise is 
significantly reduced, while the system relaxation dynam- 
ics, on the picosecond timescale, is almost unaffected by 
the steady-state thermal condition of the substrate. This 
behavior is directly connected to the rate equation gov- 
erning the temperature relaxation of the system excited 
at t=0 by the absorbed power per unit volume P(t): 

- P(t) - fcAT(t) (A2) 

being AT=T-T the temperature difference from the 
equilibrium temperature T , C the heat capacity of the 
system and k the power density dissipated per unit K 
(neglecting the temperature dependence of k). The solu- 
tion of Eq. IA2I is given by: 



i r + °° P -t'/ T 

AT(t) = - / Pit- t')- H (t')dt' + AT(0)e-'/ r iJ(i) (A3) 

k / m t 



where r—Cjk is the time constant of the system, AT(0) is the temperature difference at t—0 and H(t) is the step 
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function. 

In the limit of r— >0, being r the time constant of the 
system, the temperature variation is given by 



AT(t) 



P(t) 



(A4) 



where k is the dissipated power density normalized to the 
temperature. 

On the contrary, when r — > oo (average heating), the 
temperature variation is 

AT(t) ~ [ P(t')dt' + AT(0)H(t) (A5) 
kT Jo 

being H(t) the step function. In this case, AT(t) is pro- 
portional to the integral of the absorbed power, super- 
imposed to a constant background. In particular, if the 
excitation has a spectral component at a frequency lu, i.e. 
P(t)=P cosut, Eq. lA"5l yields to 



AT(t) 



1 sinujt 



AT(0)H(t) 



(A6) 



Therefore, if the period of the modulation of the absorbed 
power is smaller than the time necessary to the system 
to dissipate the absorbed energy, the temperature modu- 
lation at the frequency lo is damped by a factor 1 /u> and 
7r/2-out-of-phase, with respect to the excitation. 

The timescale of the relaxation dynamics of the excited 
system ranges up to the nanosecond, whereas the increase 
of the steady-state temperature of the substrate, induced 
by the modulated pump intensity P(t), lies in the mil- 
lisecond range, as discussed in Sec. IIVI Acquiring the 
100 kHz in- phase and the 7r/2-out-of-phase components 
of the probe signal, we decouple the temperature varia- 
tion in the time-scale of interest from the average heating 
of the substrate. On the contrary, standard modulation 
of the pump beam at 0.1-1 kHz (T^l-10 ms), would re- 
sult in a high in-phase background in the probe signal, 
drastically decreasing the S/N ratio. 



of the incident beam. Upon excitation with the pump 
pulse, the time-dependent variation of disk radius Sa(t) 
induces a relative variation of the reflected signal: 



51, 



refl 



(R Py R Si )2va ^ Q 28 Sa 



Irefl RSiD 2 + (RPy ~ Rsi)™ 2 

(B2) 

considering that, for our samples, D=Aa. 

Instead, the first-order diffracted electric field is deter- 
mined by the form factor (/i) of a flat disk [32| : 



J\= I e~ ir ^ds = / / re- lTUc ° s "drd6 = %yJ 1 {y) 
o Jo " 



(B3) 

being S the dot surface, r the position vector, G the re- 
ciprocal lattice vector, r the radial coordinate (parallel 
to the surface), 6 the angle between r and G and y=Ga. 
The intensity of the diffracted beam is given by the sum 
of the intensity diffracted by the dot array and, follow- 
ing Babinet's principle, of the intensity diffracted by the 
negative dot array of silicon: 

47T 2 

ho = -^{Rs*y 2 Ji{y? + R Py y 2 My?) (B4) 
With the help of the Bcsscl function relation: 

yJ[(v) = vMv)- My) (B5) 

we obtain the relative variation of the diffracted signal: 

l^ 5a = 2G ^ iRsi+Rpy)6ac , 2 . 5 5 Jt (B6) 
Iid da Mv) a 

In the diffraction configuration, the S/N ration is thus 
increased by a factor 2.5/0.28~9. Due to the fact that 
the relative intensity variation of the diffracted beam is 
of the order of 10 -5 , the increase of a factor 9 in the S/N 
ratio is necessary to obtain reliable measurements. 



APPENDIX B: DIFFRACTED INTENSITY 
VARIATION 

In this experiment the signal to noise ratio can be sig- 
nificantly improved by measuring the variation of the in- 
tensity of the first-order diffracted beam, instead of the 
reflected beam. 

The intensity of the reflected I re fi beam from a single 
unit cell of the array is given by the sum of the intensity 
reflected from the dot metal surface and of the intensity 
reflected from the complementary silicon surface, i.e. 



Irefl — 



Rsi(D 2 ~ to 2 ) + RpyTra 1 
LP 



(Bl) 



where _R.si~0.34 and Rp y ~0.62 are the silicon and 
permalloy reflectivities at normal incidence, D is the ar- 
ray period, a is the dot radius and Ii nc is the intensity 



APPENDIX C: SPACE-TIME PROFILE OF THE 
ABSORBED POWER 

We report the detailed expressions for the power ab- 
sorbed per unit volume, the boundary conditions and the 
numerical values for the parameters used in the simula- 
tions for the thermal problem addressed in section [TV] 

The space-time profile for the power per unit volume 
absorbed by the Si sample from the continous laser used 
in eq. [T] is 

P CO n(t) = 1 - exp(t/r) 
(1 — Rg-) 

Pcon(z,r) = Jo.con : exp(z/A Si )J?(r - r) (CI) 

As; 

where Io >co n = 1.68xl0 8 Wm" 2 is the equivalent average 
continuous laser power per unit surface, r=100 fs the 
laser switch-on rise-time, Asi=6.37x 10~ 6 and i?si=0.332 
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the penetration depth and normal incidence reflection co- 
efficient for Si at a wavelength of 800 nm respectively, 
calculated using the values of nsi and reported in 
Tab. A hard-edge laser beam spot size of radius ro 
= 27.5 /im is assumed, this is formally accounted for by 



TABLE III: Constants used in solving the thermal problem 
addressed in section in section IW1 



iPy 


1.065xl0 3 


(Jm" 3 K 2 ) 


[33] 


c e 


3.248xl0 5 


(Jrn^K- 1 ) 


7p h x305 K 


Cpy 


2.2xl0 6 


(Jm^K- 1 ) 


[3J 


Csi 


1.5xl0 6 


(Jm^K- 1 ) 


[35] 




60 


(Wm-'K- 1 ) 


[36] 




13-20 


(Wm-'K" 1 ) 


[3JJ 


ksi 


130 


(Wm^KT 1 ) 


[38] 


G 


8xl0 17 


(Wm^K" 1 ) 


[10] 


R th 


1Q -8 


(W _1 m 2 K) 


[13], M 


n Py 


2.25 




[34] 


k Py 


3.6 




[3J 


nsi 


3.72 




[3J 


ksi 


0.01 




[34J 


h 


5.6 


(Wm^K- 1 ) 


[39] 
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